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The  purpose  of  this  study  was  to  dsvslop  a  Method  of 
determining  a  sequence  of  optimal  orbital  evasive  and 
rendezvous  maneuvers  using  continuous,  low  thrust  electric 
propulsion.  1  began  ay  research  with  the  intent  of 
developing  an  algorithm  which  could  be  implemented  on  the 
tactical  level  to  analyze  the  optimal  avoidance  trajectory 
for  a  spacecraft  under  attack.  1  made  substantial  progress 
towards  this  goal,  and  I  believe  this  work  forms  the  basis 
of  a  practical  solution  method. 

Throughout  this  project,  the  unpredictable  behavior  of 
the  Lagrange  costates  posed  a  major  difficulty  in  that 
costate  trajectories  greatly  affect  the  convergence  of  a 
numerical  solution.  Por  this  reason,  results  Include  these 
trajectories  and  their  quadratic  curveflts  with  time  in  the 
hope  that  future  researchers  may  find  them  useful  in 
obtaining  a  converged  solution. 

Upon  completion,  I  would  like  to  add  that  this  research 
was  not  the  product  of  ay  efforts  alone.  I  owe  a  debt  of 
gratitude  to  ay  thesis  advisor,  Lt  Col  Joseph  i.  Widhala, 
for  the  endless  ideas  he  contributed  and  the  direction  he 

provided  me  over  the  course  of  ay  efforts.  Finally,  I  give  - - 

on  For 

ay  greatest  thanks  to  Keith  for  his  continuous  support  and  4&I 
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understanding  throughout  this  endeavor.  >ced 
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Aba tract 

Optimal  in-plane  orbital  maneuvers  enabling  a  target 
spacecraft  in  geosynchronous  orbit  to  avoid  interception  by 
a  threat  in  direct  coasting  ascent  are  examined.  Given 
constant  magnitude,  low  thrust  electric  propulsion,  the 
thrust  vector  control  history  is  determined  to  propel  the 
spacecraft  to  either  a  maximum  or  minimum  orbit  radius  in  a 
specified  time  of  flight.  The  evasive  maneuver  is  studied 
from  two  aspects.  First,  an  orthogonality  constraint  is 
applied  at  the  final  time  of  the  maneuver,  such  that  the 
target  spacecraft’s  position  vector  is  orthogonal  to  its 
velocity  vector.  This  requirement  is  then  removed  to  allow 
the  spacecraft  to  travel  without  influence  of  a  terminal 
constraint.  Results  demonstrate  the  adverse  shaping  effect 
the  orthogonality  constraint  applies  over  the  entire  evasive 
maneuver . 

Upon  completion  of  the  evasive  maneuver,  the  spacecraft 
returns  in  minimum  time  to  its  nominal  orbit,  in  a  position 
as  if  no  evasive  thrusting  had  occurred.  Results  show  the 
time  required  to  rendezvous  with  the  nominal  orbit  any  be 
significantly  greater  than  the  time  of  flight  of  the  evasive 
maneuver . 

Optimal  control  is  applied  to  formulate  the  evasive  and 
rendezvous  maneuvers  using  Euler -Lagrange  theory  and  the 


calculus  of  variations.  The  behavior  of  the  Lagrange 
costates  is  examined  in  terms  of  their  effect  upon  solution 
convergence.  The  resulting  two  point  boundary  value  problei 
is  solved  by  numerical  methods  to  yield  the  optimal 
spaceflight  trajectory. 
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OPTIMAL  THRUST  VECTOR  COVTROL  OP  COPLAEAR 
ORBITAL  EVASIVE  MAMEUVERS 

I .  Introduction 

In  post  ytara ,  numerous  ftudlaa  have  dealt  with  tho 
aubjoct  of  optimal  aanauvarlng  in  space ,  mainly  pertaining 
to  orbital  tranafar,  rendezvous  and  aatollita  atation 
hooping.  Until  recently ,  however,  few  authora  have 
addreeeed  the  problea  of  avoiding  an  interception  in  apace 
through  evaaive  aaneuvering.  At  thia  tiae,  anti-aatellite 
vehiclea  are  capable  of  intercepting  low  altitude  targeta 
with  no  warning  or  up  to  three  and  one  half  houra  warning. 
High  altitude  targeta  in  geoaynchronoua  orbit,  however,  aay 
receive  up  to  alx  houra  warning  of  an  Interceptor ’a  approach 
in  direct  coaating  aacent.  In  view  of  tbeae  current  and 
eaerglng  interceptor  capabilitiea ,  the  intercept  avoidance 
problew  la  becoming  increasingly  important  to  spacecraft 
adsslon  planning  in  relation  to  anti-aatellite,  apace 
reconnaissance  and  collision  avoidance  objectives. 

Much  of  the  current  literature  la  concerned  with 
optimal  maneuvering  using  impulsive  or  continuous  thrust. 
Burk  (9)  discusses  the  evaaive  problem  in  terms  of  minimum 
impulse  orbital  maneuvers.  Eel ley.  Cliff  and 
Lutxe  (8:14-17)  present  a  method  for  optimisation  of  attaak 
and  evasion  sequences  based  on  Clohessy-Wlltshire  models. 


Both  Lawdm  (0)  and  Marec  (10)  develop  thoory  directly 
applicable  to  optimal  spaceflight  trajectories ,  aithar  in 
taraa  of  intarcapt  or  avaaiva  maneuvers.  Diffarantial  gaaa 
thaory  oonaidara  thaaa  two  approachaa  simul tanaoualy  in  a 
two  apacaoraft  puraui t-evasion  acbaaa  (7) . 

Thia  work  followa  on  that  of  Praiaaingar  (13),  who 
davalopad  a  Modal  determining  orbital  avaaiva  aanauvara 
uaing  varioua  thruat  aagnitudaa.  Tha  acanario  praaantad 
harain  la  ona  side  of  tha  puraui t-avaaion  game,  with  tha 
additional  conaidaratlon  of  raturn  of  a  apacacraft  to  ita 
nominal  orbit.  A  targat  apacacraft  in  gaoaynchronoua  orbit 
ia  conaidarad  to  ba  tinder  intarcapt  by  a  non-maneuvering 
threat  in  diract  coaating  aacant.  Tha  point  of  intarcapt  ia 
apacifiad  on  tha  nominal  orbit  aithar  two,  four  or  aix  hour a 
ahead  of  the  target  apacacraft  upon  initiation  of  tha 
avaaiva  maneuver  at  a  time  t  .  Thua,  with  no  maneuvering  by 

O 

tha  targat  apacacraft,  tha  two  crafta  would  collide  at  tha 
nominal  intarcapt  time  t&.  To  avoid  any  further  af facta  of 
tha  interceptor ’ a  aalf -detonation ,  tha  targat  apacacraft 
muat  achieve  a  maximum  in-plane  change  in  orbit  radiua  auch 
that  it  will  ba  outaida  an  explosion's  lethal  radiua  at  tha 
intarcapt  time.  To  accompliah  thia  change  of  radiua,  which 
may  raault  in  aithar  a  higher  or  a  lower  altitude  trajectory 
than  tha  nominal,  tha  targat  apacecraft’a  orbit  radiua  ia 


•mm 


•  i 


maximised  or  minimised  by  applying  two  strategies  of  optimal 
control.  First,  the  targat  spacecraft  maneuvers  subject  to 


a  constraint  specifying  its  terminal  velocity  vector  as 
orthogonal  to  its  terminal  radius  vector.  The  other 
strategy  allows  the  spacecraft  to  travel  with  no  terminal 
constraint.  These  two  strategies  way  be  compared  in  terms 
of  the  effect  a  terminal  constraint  applies  over  an  evasive 
trajectory,  as  well  as  its  effect  on  the  spacecraft’s 
ability  to  return  to  its  nominal  orbit. 

The  target  spacecraft  executes  the  evasive  maneuver 
with  a  specified  amount  of  fuel,  such  that  the  maneuver 
terminates  when  the  threat  reaches  the  nominal  intercept 
point.  Considering  the  spacecraft  rocket  propellant  to  be 
consumed  at  a  small,  constant  rate  by  a  continuous  low 
thrust  propulsion  system,  the  solution  of  the  problem  is  to 
control  optimally  the  thrust  direction  with  time.  The 
position  and  velocity  of  the  target  spacecraft  are  measured 
in  a  non-rotating  x-y  coordinate  system,  oriented  such  that 
the  positive  x-axis  passes  through  the  point  where  the 
evasive  maneuver  Initiates.  The  thrust  vector  control  angle 
is  measured  relative  to  the  ♦x-axis,  as  shown  in  Figure  1. 


Immediately  upon  completion  of  the  evasive  maneuver  at 
tt,  the  spacecraft  returns  to  its  nominal  orbit.  The  return 
maneuver  is  performed  with  a  minimum  amount  of  fuel,  such 
that  the  problem  solution  is  time  optinsl .  The  evasive  and 
rendesvous  maneuvers  occur  in  sequence  over  the  time 
interval  t@  to  tf  and  remain  in  the  plane  of  the  nominal 
orbit,  as  depicted  in  Figure  2. 


For  Meh  atntuvtr  naalnad,  tho  application  of 
Kuler- Lagrange  tbaory  (4:47-80)  and  tho  calculus  of 
variations  results  in  a  two  point  boundary  value  problem 
describing  necessary  conditions  an  optimal  solution  must 
satisfy.  Tbe  differential  equations  associated  wltb  each 
optimal  control  problem  are  nonlinear,  and  tbus  must  be 
solved  by  numerical  methods. 

From  another  perspective,  the  evasive  maneuver  may  be 
viewed  in  the  context  of  intercepting  a  specified  radius  in 
minimum  time.  Two  strategies  are  again  implemented,  such 
that  the  terminal  orthogonality  constraint  is  first  applied 
and  then  removed.  This  viewpoint  is  considered  to  ensure 
the  numerical  solutions  aquired  for  the  evasive  maneuvers 
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Optical  control  problems  ouch  an  those  discussed  in 
in  thin  study  art  dovolopod  using  tho  calculus  of 
variations.  Bryson  and  Bo  (4:47-80)  discuss  tbs  following 
method  of  applied  optical  control  as  tbs  Kulsr-Lagrangs 
equations.  Tbis  theory  completely  defines  a  two  point 
boundary  value  problem  to  be  solved  by  analytical  or 
numerical  means,  with  either  a  fixed  or  free  terminal 
time  tf. 

Tbe  system  under  consideration  is  described  by  a  set  of 
nonlinear  differential  equations 

x  «  f[x(t> ,t]  (1) 

where  the  states  are  specified  at  an  Initial  time  t^. 

A  performance  index,  composed  of  functions  to  be 
maximised  or  minimised,  is  written  in  the  scalar  form 


*(x(tf).tf]  ♦  j|f  Jf[x(t)  ,fHt)  ,t]dt 


A  set  of  functions  of  the  terminal  state  aonstrains  the 


system  at  the  final  time: 


Adjoining  these  constraint*  and  tha  system  state  equations 
(1)  to  the  perforaance  index  (2)  by  a  set  of  constant 
multipliers  v  and  costate  multipliers  X  yields 


J  «  g>[x(tf).tf]  ♦  uTjg[x(tf)  ,tfl  ♦  J^f  (JfCx(t)  ,f?(t)  ,tl 


♦  X  (f  -  x))dt 


(4) 


A  dimensionless  scalar  function  X,  the  Hamiltonian,  is 
defined  as 


*[x(t)  ,X(t)  ,tl  =  Jffx(t)  .ftlt)  ,t] 

♦  ^T<t)f  Cx(t)  ,/?(t)  ,t]  (5) 

Substituting  for  X  in  Eq  (4)  and  taking  the 
differential  to  allow  for  infinitesimal  changes  in  the 
terminal  time  t, ,  considering  t  fixed,  yields 

I  O 


Integrating  the  last  term  on  the  right  hand  side  of 
*q  (0)  by  parts  and  substituting  for  the  variation  of  x  in 
terms  of  the  differential  dx 
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c\  1,U 


to  combine  to rms  in  dx  gives 


{[  %  ♦  x*SE  *  *  ♦  ]«,  *  [  %  '*'%  -  iT 

♦  <hTsi>  l'  *  J/  {[  If  *  iT]*S  *  If  *«J«‘  <8 


Choosing  the  Multiplier  functions  X(t)  to  cause  the 
coefficients  of  <5x(t)  ,  dx(tf)  end  dtf .  if  the  final  t ime  is 
undetermined ,  to  vanish  gives  a  set  of  conditions  necessary 
for  the  existence  of  a  stationary  value  of  J.  These 
necessary  conditions  include  an  optimality  condition 
determining  the  control  ft{t) 


as  well  as  a  set  of  Lagrange  costate  equations 


vith  boundary  conditions  given  by 


/«.« .  r £ .  U 


•!  * 


An  additional  transversality  condition  applies  if  the 
terminal  tine  tf  is  unspecified: 


r  *2 

[  at 


+ 


0 


(12) 


On  the  optimal  trajectory,  Eq  (9)  shows  that  the 
optimal  control  ft( t)  maximizes  or  minimizes  the  Hamiltonian. 
It  then  follows  that  to  maximize  X,  and  thus  maximize  J,  the 
following  statement  known  as  the  Legendre-Clebsch  condition 
in  the  calculus  of  variations  must  be  enforced: 


aft2 


<  o 


(13) 


To  minimize  0f,  and  thus  J,  however,  the  Legendre-Clebsch 
condition  must  take  the  form 


o2x 

aft2 


>  0 


(14) 


Thus,  the  two  point  boundary  value  problem  is 
completely  described  by  the  state  equations  (1)  with  the 
initial  states  specified,  and  the  costate  equations  (10) 
with  boundary  conditions  given  as  in  Eq  (11).  The  control 
is  determined  by  the  optimality  condition  (9) ,  subject  to 
the  Legendre-Clebsch  condition  of  Eq  (13)  or  (14).  The 
remaining  boundary  conditions  are  recovered  from  the 
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A  target  spacecraft  is  considered  to  be  threatened  by 
a  coasting  craft  in  a  dihect-ascent  collision  course  with 
its  nominal  geosynchronous  orbit.  To  avoid  intercept  by  the 
threat  craft  and  any  possible  effects  of  its 
self -detonation,  the  target  craft  must  maneuver  in  an 
optimal  manner  to  achieve  a  maximum  change  in  its  orbit 
radius.  This  change  in  the  spacecraft’s  trajectory  may  be 
directed  to  the  exterior  of  its  nominal  course,  or  to  the 
interior  providing  it  does  not  pass  within  a  lethal  radius 
centered  on  the  ascending  interceptor.  In  the  context  of 
optimal  control,  these  two  types  of  trajectory  changes  are 
effected  by  maximizing  and  minimizing  the  entire  radius  of 
the  target  spacecraft's  orbit. 

The  evasive  maneuver  initiates  at  time  t^,  prior  to  the 
nominal  intercept  time  t&  by  a  specified  time  of  flight. 
Continuous,  low  thrust  electric  propulsion  is  employed  over 
the  entire  naneuver ,  wherein  control  is  provided  by  a  thrust 
vector  free  to  gimbal  in  any  desired  direction.  Both  the 
maximum  and  the  minimum  radius  maneuvers  are  accomplished  in 
two  ways.  First,  the  target  spacecraft’s  orbit  radius  is 
maximized  under  the  influence  of  a  terminal  orthogonality 
constraint,  such  that  the  spacecraft's  velocity  vector  is 
orthogonal  to  its  position  vector  at  the  final  time  of  the 
maneuver,  t  .  The  constraint  is  then  removed  to  allow  the 
spacecraft  to  travel  unhindered  by  any  shaping  effect  the 


terminal  condition  might  have  on  the  maneuver.  These  two 


classes  of  evasion  maneuvers  may  be  compared  in  terms  of 
total  fuel  required  to  achieve  a  maximal  radius  change  at  t^ 
followed  by  a  rendezvous  with  the  nominal  orbit  at  time  tf. 

Equations  of  Motion 

The  spacecraft  travels  in  a  plane  surrounding  the 
earth,  such  that  its  position  is  described  in  non-rotating, 
earth-centered  coordinates  x  and  y.  The  spacecraft  velocity 
components  are  given  by  u  and  v  along  the  x  and  y  axes, 
respectively.  The  two  body  equation  of  motion  for  a 
spacecraft  operating  in  an  inverse  square  gravitational 
field  without  influence  of  external  forces  is  developed  by 
Bate  (3:14)  and  Baker  (1:58).  Veglecting  perturbative 
effects  and  taking  the  rocket  thrust  force  into 
consideration,  the  equations  of  motion  may  then  be  written 
in  first  order  form  as  given  by  Preissenger  (13:7) : 


x  *  u  (15) 

y  =  v  (16) 

u  *  -  ♦  - — — —  co mft  (17) 

r  m  -  mt 

O 

v  *  ♦  - T  ;  -  sin/3  (18) 

r  m  -  mt 

O 

where  ft  is  the  thrust  vector  control  angle  measured  with 


respect  to  the  *x-axis.  Assuming  the  orientation  of  the 


x-axia  ia  described  as  passing  through  the  point  where  the 


$(x(t&)  ,ttJ  «  r(tt)  *  (x*  ♦  y2)1/2|t  (20) 

*  t 

Adjoining  the  constraint  of  Bq  (10)  to  Bq  (20)  yialda  a 
performance  index  posed  as  in  Bq  (4) : 

J  *  r(t4)  ♦  v  (xu  ♦  yv)  I  (21) 

*  > 

ehere  the  multiplier  u  is  to  be  determined  by  the  boundary 
conditions  on  the  Lagrange  costates.  The  problem  is  now  to 
determine  the  thrust  vector  control  angle  ft  to  maximize  this 
performance  index. 

Following  the  Lagrange  method  described  by  Bryson  and 
Ho  (4) .  the  system  Hamiltonian  is  constructed  as  given  by 
Bq  (5)  : 


X 


X  u  ♦  X  v 

x  y 


- ; —  com  ft 

m  -  mt 


♦ 


(22) 


The  thrust  control  angle  producing  a  stationary  value 


of  the  performance  index  Bq  (21)  is  given  by  the  optimality 
condition  of  Bq  (0),  resulting  in  a  bilinear  tangent  control 


To  determine  the  Lagrange  multipliers  X,  the  coetate 


equation*  are  developed  from  Sq  (10)  as 


X 

V 


-X 


y 


(28] 
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(3o: 

(31 


leceasary  boundary  conditions  on  the  coatates  at  the  final 
time  of  the  evasive  auuieuver  are  determined  from  Sq  (11): 


X  (t  ) 

X  1 


*  ♦  vu 
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X  (t  )  1  -  ♦  vv 
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X  (t)  = 

u  t 
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X  (t  ) 

V  1 


(32: 
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(34 

(35 


Combining  these  boundary  conditions  with  the  orthogonality 
constraint  function  of  Sq  (19)  as  moll  as  the  specified 
Initial  state  values  to  eliminate  the  multiplier  v  yields  a 
set  of  boundary  conditions  implemented  to  solve  the 


resulting  two  point  boundary  value  problem: 


x(t  )  Spocifiod 

~  O 

xu  ♦  yv  *  0  |t 


X.  x  -  X  y 

V  u 


X  X 

X 


X  y  - 
y 


(30) 

(37) 

(38) 

(30) 


Theao  boundary  conditions  (30)  through  (30)  art  not  unique 
in  form,  but  they  must  be  independent  at  time  t  . 

The  two  point  boundary  value  problem  is  now  described 
by  the  four  state  equations  (IS)  through  (18)  with  the 
appropriate  definitions  of  the  control  angle  applied, 

Eqs  (25)  and  (20) ,  the  four  costate  equations  (28)  through 
(31),  four  specified  initial  states,  and  the  above  four 
final  boundary  conditions  on  the  costates. 

Mlniaua  Radius  Maneuver .  Depending  upon  the  approach 
path  of  the  attacker,  an  avoidance  saneuver  to  the  interior 
of  the  nominal  spacecraft  orbit  way  be  practical.  Such  a 
trajectory  is  best  Modeled  as  a  Minimum  radius  maneuver  or, 
alternatively ,  a  maximum  negative  radius  maneuver.  The 
development  of  the  two  point  boundary  value  problem  defining 
the  maneuver  follows  that  of  the  maximum  orbit  radius 


■tntuvtr ,  with  wnly  a  ai|n  ohtnga  in  the  boundary 
conditions . 

To  maximise  tha  negative  of  the  spacecraft  orbit 
radius,  subject  to  the  orthogonality  constraint  on  the 
spacecraft  velocity  and  position  vectors  at  the  final  tiwa 
the  performance  index  is  written  as 

J  =  -rCt^)  ♦  w(xu  ♦  yv)  I  (4 

•  i 

Since  the  performance  index  is  again  to  be  maximized, 
the  Hamiltonian  Kq  (22) ,  control  law  Kq  (23) . 
Legendre-Clebsch  condition  Hq  (24) ,  control  angle 
definitions  Eqs  (29)  and  (20).  and  costate  equations  (28) 
through  (31)  determined  previously  apply.  However,  the 
necessary  boundary  conditions  on  the  costate  equations 
are  now 
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Kliminatlng  the  multiplier  v  and  considering  the 


orthogonality  constraint  Kq  (10)  yields  tho  following  sot  of 
boundary  conditions  to  apply  to  tho  two  point  boundary  valuo 
problow: 


g(to>  Spocifiod 
xu  ♦  yv  ■  0 


X  x  -  X  y 

V  u 


k 


X  x  ♦ - X  u  *  0 

x  r  u 


Xy  +  ^-  Xw^O 
y  r  v 


(43) 

(40) 

(47) 

(48) 


Ago  in,  thoso  boundary  conditions  aro  indopondont  at  tiwo  t  . 

Maneuver  with  no  Torwinal  Constraint 

In  sows  situations,  tho  spacocraft  way  rocoivo  little 
warning  of  an  attacker’s  approach.  It  way  then  bo  necessary 
for  tho  spacocraft  to  travel  as  far  away  frow  tho  throat  as 
possible  in  a  specified  tiwo  without  regard  for  tho  final 
oriontation  of  its  velocity  vector.  In  effect,  this  allows 
the  target  spacecraft  to  achieve  the  saws  radial  distance 
change  as  achieved  with  a  torwinal  constraint  in  a  wuch 
shorter  tiwe  period. 

Masiwuw  Radius  Maneuver .  The  requirewent  of  achieving 
waxlwuw  radius  is  still  expressed  by  the  function  of 


Bq  (20).  In  this  case,  however,  the  perforwance  index  to  be 


The  definitions  of  ths  control  angle  ft  mist  then  be 
converted  at  the  final  tin  from  Eqs  (25)  and  (26)  to 


cos/?(  t^) 
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The  Legendre-Clebach  condition  is  then  satisfied  at  the 
final  tine  as 


* - - -  (X  2  ♦  X  2)1/2  I  (57) 

1  m  -  it  *  y  It 

o  1 

The  two  point  boundary  value  problem  is  now  given 
completely  by  the  state  equations  (15)  through  (18)  with  the 
appropriate  control  definitions  applied  over  the  maneuver 
and  at  the  final  time,  the  costate  equations  (28)  through 
(31) ,  the  boundary  conditions  on  the  costates  at  the  final 
time,  Eqs  (50)  through  (53),  and  the  specified  states  at  the 


Initial  ti 


IV.  Mini mum  Time  Bendezvoua  Maneuvers 


Once  the  target  spacecraft  has  completed  an  evasive 
manevuer  at  time  t  ,  the  concern  is  to  return  the  spacecraft 
to  its  nominal  orbit  in  a  time  optimal  manner.  Ideally,  the 
point  of  rendezvous  is  the  nominal  geosynchronous  position 
of  the  target  spacecraft  had  it  never  performed  the  evasive 
maneuver . 

Equations  of  Motion 

The  state  equations  (15)  through  (18)  again  govern  the 
spacecraft  motion  with  time.  However,  in  order  to 
numerically  determine  the  final  optimal  time  of  the 
rendezvous  maneuver  over  a  fixed  time  interval,  a  change  of 
variables  is  perforsmd  such  that 

t  =  ct  (63) 

where  t  is  defined  over  the  interval  T  =  0  to  T  =  1 
c  is  then  a  time  parameter  state  with  a  constant  value  equal 
to  tf,  the  final  time  of  the  return  maneuver.  Applying  the 
transformation 


d  (  )  _  d(  ) 
dT  c  dt 


(04) 


to  the  state  equations  yields  the  parameterized  state 


equations : 


where  the  prime  symbol  indicates  differentiation  with 


respect  to  t.  The  time  parameter  adds  an  additional  state 
equation : 


c'  =  0  (09 

Assuming  the  rendezvous  maneuver  initiates  immediately 

following  completion  of  an  evasion,  the  initial  conditions 

given  to  the  state  variables  x(t  )  are  the  terminal  state 

~  1 

values  determined  at  the  final  time  of  the  avoidance 
maneuver . 

Minimum  Time  Maneuver 

The  rendezvous  is  performed  in  minimum  time,  such  that 
the  performance  index  is  written  from  Eq  (2)  as 

J  = 


I' 


ldt 


(70 


Applying  Eq  (04)  to  the  above  yields  the  transformed 


V  =  V  COS 

NOM  O 


i  L  ‘  *  -  v  J  j 


In  thtaa  equations,  r  ia  the  radius  of  the  nominal 

O 

geosynchronous  orbit,  is  the  nominal  circular  speed,  and 
(t^  -  tQ)  is  the  time  of  flight  of  the  original  evasive 
maneuver . 

Again  parameterizing  time  as  done  for  the  performance 
index,  the  nominal  orbit  equations  becos» 
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The  constraints  of  Eqs  (72)  through  (75)  are  then  evaluated  at 
the  final  tine  r  =  1 

Adjoining  the  constraint  equations  to  the  performance 
index  of  Eq  (71)  leads  to  an  expression  to  be  minimized: 


J 


=  I  V  (  x  -  x _ )  +  V  (  y  -  y  )  ♦  V  (  u  -  u  ) 

|  X  NOM  y  9  9  NON  U  NOW 

T  =1  Jo 


♦  V  (  v  -  V  ) 
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cdT  (84) 


The  system  Hamililtonian  is  now  assembled  as  in  Eq  (5) 


flf=c-(l  +  ^u  +  \v  ♦  \ 
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ZJ^~  *  - sinf?  |  ^  (85) 
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The  performance  index  J  is  again  given  a  stationary 
value  subject  to  the  optimality  condition  of  Eq  (8) .  As 
pertained  to  the  evasive  maneuvers,  the  bilinear  tangent 
control  law  applies,  but  the  thrust  control  angle  definition 
changes  such  that 


ft  * 


(86) 


This  change  is  attributed  to  the  Legendre-Clebscb  condition, 
which,  from  Eq  (14),  for  a  minimization  problem  requires 
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nn 
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To  satisfy  this  condition,  tho  thrust  control  angle  Bust  bo 
appropriately  defined  in  the  state  equations  (85)  through 
(88)  as: 


ft 


cos/?  =  - 


.  2  .  2 

X.  ♦  X 

U  V 


sin/?  =  - 


X  *  v  2 

X  ♦  X 

U  V 


thus  yielding  the  control  law  given  in  Eq  (86)  and 
satisfying  the  Legendre-Clebsch  condition  as 
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X  2  ♦  X  * 

U  V 


The  costate  equations  are  transforwed  fros  those  of 
Eqs  (28)  through  (31)  as 
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Th«  boundary  conditions  on  tha  coatatas  at  tha  final 


ar(T=l)  ♦  X  y  -  X  x  ♦  X  v 

[  x'NOU  y  NOH  u  NOW 


-Xu 
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=  0  (100) 


■her*  the  Hamiltonian  X  is  given  by  Eq  (85)  and  the  nominal 
orbit  states  are  given  by  Eqs  (76)  through  (70) .  It  should 
again  be  noted  that  the  final  t i me  tf  corresponds  to  a  value 
of  r  a  1  Additionally,  the  constant  Lagrange 

multipliers  v  mere  eliminated  from  Eq  (100)  by  applying 
boundary  conditions  on  the  costates  given  by  Eq  (95) . 

Thus,  the  tmo  point  boundary  value  problem  is 
completely  defined  by  the  state  equations  (65)  through  (68), 
the  costate  equations  (91)  through  (94) ,  four  boundary 
conditions  specifying  the  initial  state  at  tjt  four  boundary 
conditions  on  the  final  state  at  tf,  and  the  transversal ity 
condition  of  Eq  (99) . 


V.  verification  of  Fixed  Time  Solution* 


Since  the  maximum  radius  evasive  maneuvers  yield 
optimal  trajectories  over  a  specified  time,  the  two  point 
boundary  value  problem  describing  a  minimum  time  intercept 
of  a  specified  radius  should  ideally  generate  an  identical 
solution.  The  formulation  of  the  time  optimal  intercept  is 
similar  to  that  of  the  rendezvous  problem,  and  it  provides  a 
convenient  swans  of  verifying  the  evasive  maneuver  control 
history . 


Maneuver  with  Terminal  Constraint 


The  requirement  of  intercepting  a  specified  orbit 
radius  at  the  final  time  of  the  maneuver  t^  is  expressed 
through  the  constraint  function 


Vjt^xU^  J  =  x(tt)2  ♦  y(tt)2  -  r(tt)2 


(101) 


where  r(t^)  is  the  maximized  orbital  radius  achieved  in  the 
evasive  ausneuver.  Similarly,  the  orthogonality  requirement 
is  defined  through  the  constraint 


C t4,x(t  )  ]  *  xu  ♦  yv 


(102) 


By  adjoining  the  constraint  equations  (101)  and  (102)  to  the 
minimum  time  performance  index  ,  Eq  (71) ,  with  a  pair  of 


multipliers  u ,  the  two  point  boundary  value  problem  la 
davalopad  to  minimize  tha  new  perforaanca  indax 

J  *  i^fx2  ♦  y2  -  r2)  |  v^(xu  ♦  yv)  J  ♦  J  cdT  (103) 

T-i  T-i  0 

where  time  has  baan  paraaetarlzed  through  E q  (04) . 

Since  the  formulation  of  tha  intercept  problem  pertains 
to  minimum  time,  tha  state  aquations  (05)  through  (08), 
Hamiltonian  Eq  (85) ,  control  law  Eq  (80) ,  Legendre-Clebsch 
condition  Eq  (87),  control  angle  definitions  Eqs  (88)  and 
(80) ,  and  costate  equations  (01)  through  (04)  of  the 
rendezvous  problem  apply.  However,  the  necessary  boundary 
conditions  on  the  costates  at  the  final  time  t±,  are  changed 
to  give,  from  Eq  (11) 


X  (t*1)  ■  2v  x  ♦  v  u  (104) 

x  12 

X  (r*l)  *  2u  y  ♦  v  v  (105) 

y  1  2 

X  (T-l)  *  U  x  (108) 

u  S 

X  (T-l)  -  V  y  (107) 

v  2 


Since  the  final  time  of  the  maneuver  is  unspecif led ,  an 
additional  boundary  condition  is  formulated  through  the 
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tranavtrial ity  condition  of  Eq  (12)  as 

ar(r*i)  =  o  doe) 

Khar*  the  Hamiltonian  is  again  given  by  Eq  (85). 

Eliminating  the  multipliers  v  in  the  costate  boundary 
conditions  and  recovering  the  remaining  boundary  conditions 
from  the  constraint  equations  (101)  and  (102)  results  in  a 
set  of  boundary  conditions  implemented  to  solve  the  tmo 
point  boundary  value  problem.  These  conditions  are 


x(t=0)  Specified 

xu  ♦  yv  *  0  1 

■  r=i 

x  ♦  y  -  r  *  0  S 

•  T  =1 

X  y  -  X  x  *  0  I 

u  v  I 

*T»1 

x(xX  -  vX  )  ♦  y(x iX  -  xX  )  =0 

y  u  u  x 

and  the  transversal ity  condition  Eq  (108).  Mote  that  Eqs 
(111)  and  (112)  are  not  the  only  means  of  combining  the 
costate  equations  to  eliminate  the  multipliers  v .  In 
general,  the  form  of  the  boundary  conditions  is  not  unique, 
but  the  conditions  must  be  independent. 
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(109) 

(110) 

(111) 

(112) 
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The  two  point  boundary  value  problem  ia  now  completely 


i-H 


defined  by  the  state  equations  (05)  through  (6B)  ,  with 
initial  conditions  specified  at  the  starting  point  of  the 
evasive  maneuver,  the  costate  equations  (01)  through  (04), 
the  transversal ity  condition  Eq  (108),  and  the  boundary 
conditions  Eqs  (100)  through  (112).  The  control  angle 
definition  of  Eqs  (88)  and  (80)  apply  to  the  state 
equations . 


Maneuver  with  no  Terminal  Constraint 

To  verify  the  unconstrained  maximum  radius  maneuvers, 
the  orthogonality  constraint  is  removed  from  Eq  (103)  to 
give  the  performance  index 


2  2  2  I  f4 

J  ■  ^(x  ♦  y  -  r  )  |  ♦  J 


(113) 


Again,  the  Hamiltonian  of  Eq  (85),  control  law  Eq  (86), 
Legendre-Clebsch  condition  (87),  control  angle  definitions 
Eqs  (88)  and  (80)  ,  and  costate  equations  (01)  through  (04) 
of  the  rendezvous  problem  apply.  In  this  case,  the  necessary 
boundary  conditions  on  the  costates  at  the  final  time  t  ,  or 
t  s  1  ,  are  derived  from  Eq  (11)  as 


\  (t*1)  ■  2v  x 

x  1 


X.  (T=l)  *  2v  y 
y  * 


(114) 


(115) 
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X  (T=l)  =  0 
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X  (T=l)  =0  (1 


As  m a  t ound  for  the  unconstrained  maximum  radius  maneuvers 
the  boundary  conditions  on  the  X  and  X  costates  yield  a 

U  V 

singularity  in  the  control  law  of  Eq  (86)  at  the  final  time 
Applying  L’ Hospital* s  rule  and  the  costate  equations  (93) 
and  (94) ,  the  control  law  at  t^  can  be  found  as 


/?(t=  1 )  *  tan  1 

1 - 

/<>U‘ 
1  1 

*  tan  1 

■  -X  ' 

X 

-X" 

L  “  J 

T  =1 

y 

(118) 


which,  when  combined  with  the  Legendre-Clebsch  condition, 
gives  the  following  control  angle  definitions  applicable  to 
the  state  equations  at  the  final  time: 
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COS^(T=l) 
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thus  satisfying  the  Legendre-Clebsch  condition  at  t  *  1  ai 
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Th«  two  point  boundary  value  problem  is  now  completely 
defined  a a  for  the  constrained  intercept  maneuver,  with  the 
exception  of  directly  applying  the  costate  boundary 
conditions  (114)  through  (117)  and  converting  the  control 
angle  definitions  to  those  of  Sqs  (110)  and  (120)  at  the 


final  ti 


I 


VI .  Method  of  Solution 

The  formulation  of  the  two  point  boundary  value  problem 
using  optimal  control  theory  produces  a  set  of  nonlinear 
ordinary  differential  equations  which  must  be  solved  by 
numerical  methods.  Each  of  the  two  point  boundary  value 
problems  developed  previously  were  solved  directly  using  the 
IMSL  routine  DVCPR.  The  algorithm  applied  in  DVCPR  is  a 
variable  order,  variable  step  size  finite  difference  method 
with  deferred  corrections,  based  upon  Lentini  and  Pereyra's 
PASVA3  program  (12:67-88).  DVCPR  solves  two  point  boundary 
value  problems  over  a  fixed  time  interval,  thus  the 
transformation  over  the  interval  t  =  0  to  t  -  1  applied  in 
the  minimum  time  maneuvers. 

As  shown  by  Carnahan,  Luther  and  Wilkes  (6:405-506), 
no  known  numerical  methods  exist  to  guarantee  a  solution  for 
a  given  boundary  value  problem.  Bryson  and  Ho  (4:214-215) 
additionally  comment  that  small  errors  in  guessing  initial 
conditions  often  produce  extremely  volatile  state  space 
trajectories,  wherein  the  values  of  the  states  and  costates 
may  well  exceed  the  computer’s  numerical  limits.  The 
greatest  difficulty  in  producing  an  accurate  solution  to  a 
well-posed  two  point  boundary  value  problem  is  then  to 
provide  the  algorithm  with  a  trial  solution  suitably  close 
to  satisfying  the  differential  equations  and  boundary 


conditions . 


The  IMSL  routine  DVCPB  demands  input  of  an  initial  grid 
to  begin  Vewtonian  iteration  on  estimates  of  unspecified 
initial  or  terminal  conditions  in  order  to  satisfy  the 
differential  equations  and  specified  boundary  conditions. 
This  grid  consists  of  trial  state  and  costate  trajectories 
over  the  entire  time  interval  of  a  maneuver.  The  input  grid 
for  the  maximum  radius  maneuver  with  the  orthogonality 
constraint  applied  at  a  final  time  of  six  hours  was 
constructed  first,  since  output  for  a  similar  maneuver  was 
available  from  the  work  of  Preissinger  (13).  Curvefits  with 
time  over  both  the  states  and  the  costates  of  Preissinger ’ s 
solution  provided  a  suitable  grid  allowing  DVCPB  to  converge 
upon  a  solution.  The  form  in  which  the  final  boundary 
conditions  were  specified  was  also  a  consideration  since, 
by  trial  and  error,  some  forms  were  found  to  be  more 
conducive  to  convergence  than  others.  For  maximum  radius 
maneuvers  over  the  two  hour  and  four  hour  times  of  flight 
with  the  orthogonality  constraint  applied  at  the  final  time, 
the  six  hour  grid  was  scaled  over  the  appropriate  time 
interval  before  input  to  achieve  a  converged  solution. 

Sample  subroutines  called  by  DVCPB  for  these  maneuvers  are 
listed  in  appendix  A. 

Maximum  radius  maneuvers  with  no  constraint  at  the 
final  time  were  considered  next,  since  these  involved  merely 
eliminating  portions  of  the  final  boundary  conditions 
specified  in  the  previous  maneuver.  The  input  grid  for 


each  time  of  flight  was  the  output  for  the  same  time  of 
flight  with  the  tangency  constraint  applied  at  the  final 
time.  Appendix  B  contains  subroutines  called  by  DVCPB  for 
these  maneuvers . 

A  similar  procedure  was  followed  for  the  mini mum  radius 
maneuvers.  Again,  the  grid  input  to  the  algorithm  was  the 
output  grid  for  the  same  time  of  flight  from  the  maximum 
radius  maneuver  with  the  corresponding  constraints  at  the 
final  time,  an  orthogonality  constraint  or  no  constraint. 

For  the  minimum  time  intercept  maneuvers  used  to 
confirm  the  maximum  radius  results,  the  grid  was  normalized 
over  time  to  achieve  convergence.  Additionally,  due  to  the 
change  in  control  angle  definition  derived  from  the 
Legendre-Clebsch  condition,  a  change  in  sign  of  two  of  the 
costates  was  required.  Again,  boundary  conditions  were  a 
consideration  in  that  the  form  they  were  written  in 
influenced  the  ability  of  DVCFR  to  produce  a  converged 
Solution . 

Rendezvous  maneuvers,  however,  posed  a  difficulty  since 
the  states  were  shifted  greatly  from  those  found  in  the 
evasive  maneuvers.  In  this  case,  the  input  grid  consisted 
of  states  approximated  in  the  appropriate  nominal  orbit 
position,  and  costates  output  from  the  solution  of  the 
minimum  time  intercept  maneuvers.  The  rendezvous  position 
was  initially  specified  as  a  stationary  point  on  the  nominal 
orbit.  A  continuation  scheme  was  then  applied  such  that  the 


final  boundary  conditions  were  eoved  from  a  stationary 
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VII.  Results 


Results  were  calculated  for  a  nominal  2400  kilogram 

spacecraft  propelled  by  ten  30  centimeter  diameter  primary 

propulsion  electric  thrusters  normally  used  for  satellite 

stationkeeping  or  interplanetary  spaceflight.  The 

propulsion  system  delivers  1.3  Hewton  total  thrust,  with  a 

-s 

total  mass  flow  rate  of  6.0  x  10  kilogram  per  second 
(11:14-10) . 

Numerical  solutions  for  each  optimal  control  problem 
were  generated  using  the  IMSL  routine  DVCPR  on  the  ASC 
VAX-11/785  system.  All  optisal  trajectories  were  determined 
to  an  accuracy  within  a  range  of  10  11  to  10  12 ,  as 
specified  by  the  tolerance  parameter  internal  to  the  IMSL 
routine  DVCPR.  The  grid  sizes  required  by  DVCPR  for 
convergence  varied  between  200  and  300  subintervals  for  the 
evasive  maneuvers.  Rendezvous  aaneuvers ,  on  the  other  hand, 
required  grid  sizes  on  the  order  of  300  to  600  subintervals. 
The  small  tolerance  parameter  and  large  grid  sizes  give  an 
indication  of  the  quality  of  the  solution. 


Control  Histories  and  State  Trajectories 

Por  all  maneuvers  with  a  specified  terminal  time,  two, 
f our  and  six  hour  times  of  flight  were  considered.  During 
these  fixed  time  maneuvers,  the  nominal  orbit  traveled 
one-twelfth,  one-sixth  and  one-quarter  of  a  geosynchronous 
orbit.  Kach  maneuver  consumed  less  than  0.1  percent  of  the 


total  fuel  Mass  allotted  to  the  spacecraft,  such  that  the 
target  spacecraft  could  maneuver  over  the  entire  design 
lifetime  of  the  propulsion  system,  if  necessary. 

Evasive  Maneuver  with  Terminal  Constraint .  Figures  3, 

4  and  5  show  the  thrust  vector  control  histories  for  two, 
four  and  six  hour  time  of  flight  maneuvers,  respectively,  to 
a  maximum  orbit  radius  with  the  orthogonality  constraint 
applied  at  the  final  time.  The  two  and  four  hour  thrust 
angle  profiles  show  similar  control  characteristics;  the 
spacecraft  initially  thrusts  into  the  first  quadrant  to 
increase  its  radial  position,  as  expected.  Mote,  however, 
that  the  terminal  control  angle  is  in  the  third  quadrant,  in 
effect  reversing  the  thrust  direction  to  achieve 
orthogonality  at  the  final  time  of  the  maneuver.  The 
control  angle  behaves  in  a  different  manner  over  a  six  hour 
time  of  flight.  Initially,  the  spacecraft  thrusts  into  the 
second  quadrant,  anticipating  its  terminal  position  nearly  a 
quarter  orbit  away.  The  thrust  angle  gimbals  in  a  direction 
opposite  to  that  of  the  two  and  four  hour  time  of  flight 
maneuvers ,  and  is  directed  nearly  along  the  negative  y-axis 
at  the  final  time  in  order  to  achieve  orthogonality. 

For  all  three  constrained  evasive  maneuvers,  the 
control  characteristic  of  reversing  thrust  direction  to 
satisfy  the  tangency  constraint  at  the  final  time  shapes  the 
maneuver  over  the  entire  time  of  flight.  This  shaping 
adversely  affects  the  increase  in  orbit  radius  achieved  at 


w 

.  A 

vv 


A1 


1 


■  Fi* 

.'-V* 

»v 


43 


THRUST  CONTROL  ANGLE  (DEGREES) 


the  final  time,  a«  will  be  seen  in  the  discussion  of  the 


unconstrained  evasive  maneuver.  The  six  hour  evasive 
maneuver  also  shows  a  unique  feature  in  that  it  terminates 
on  the  +y-axis,  on  the  same  radius  vector  as  the  nominal 
orbit  and  thus  retaining  the  nominal  orbit  period.  In 
contrast,  the  two  and  four  hour  maneuvers  achieve 
a  greater  relative  increase  in  orbit  radius  at  the  expense 
of  orbit  drift,  an  expected  phenomenon  wherein  the  orbital 
period  of  a  spaceflight  trajectory  increases  with  increasing 
orbit  radius.  As  a  result,  at  the  terminal  time  of  the  two 
and  four  hour  evasive  maneuvers,  the  target  spacecraft's 
siotion  lags  behind  the  nominal  orbit.  The  effect  herein  is 
evidenced  in  the  slight  radial  distance  benefit  gained  by 
increasing  the  time  of  flight  from  four  to  six  hours,  as 
opposed  to  the  great  amount  gained  by  increasing  the  time  of 
flight  from  two  to  four  hours.  Table  I  summarizes  the 
distance  increases  achieved  in  the  radial  direction  for  each 
time  of  flight. 


Table  I.  Radial  Increase  for  Maximum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of  Flight  (hrs) 

Radial  Increase  (km) 

2 

14. 1625 

4 

04.4401 

e 

71.0190 

WW 


For  the  purpose  of  verifying  boundary  conditions  and  the 


orthogonality  constraint.  Tables  II  and  III  give  the 
terminal  position  and  velocity  states  of  the  target 
spacecraft  upon  completion  of  the  evasive  maneuver. 


Table  II.  Terminal  Position  States  for  Maximum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of  Flight 

(hrs)  x  (km) 

y  (km) 

2 

36600.8575 

21116.0609 

4 

21216.5804 

36600.7937 

6 

o 

o 

42312.7415 

Table  III.  Terminal  Velocity  States  for  Maximum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 

Time  of  Flight 

(hrs)  u  (km/s) 

v  (km/8) 

2 

-1.5335 

2.6580 

4 

-2.6501 

1.5362 

6 

-3.0691 

0.0 

Figures  6,  7  and  8  illustrate  the  control  histories  of 
minimum  radius  evasive  maneuvers  for  two,  four  and  six  hour 
times  of  flight,  respectively.  Each  thrust  angle  profile 
shows  the  same  characteristics  as  the  parallel  maximum 
radius  control  history,  however  the  control  angle  is  shifted 
by  180  degrees.  This  shift  results  in  nearly  identical 
changes  in  orbit  radius  at  the  terminal  tisw,  but  in  a 


radially  opposite  direction.  Table  IV  lists  the  results. 


Table  IT.  Badial  Decrease  for  Minimum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of  Flight  (hrs) 


2 


Radial  Decrease  (km) 


14.1632 
64.5499 
72.0168 


The  final  position  and  velocity  states  for  the  constrained 
minimum  radius  maneuver  are  given  in  Tables  T  and  ¥1 . 


Table  T.  Terminal  Position  States  for  Minimum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of  Flight  (hrs) 


2 


x  (km) 


y  (km) 


36562.9713 

21124.9903 

21204.5345 

36562.7162 

o 

• 

o 

42169.1057 

Table  VI .  Terminal  Velocity  States  for  Minimum  Radius 
Evasive  Maneuver  with  Terminal  Constraint 


Tims  of  Flight  (hrs)  u  (km/a) 


v  (km/s) 


Evasive  Maneuver  with  no  Terminal  Constraint 


thrust  angle  profiles  with  time  for  the  two,  four  and  six 
hour  time  of  flight  maneuvers  to  a  maximum  orbit  radius 
without  terminal  constraint  are  shown  in  Figures  9,  10  and 


11,  respectively.  Due  to  the  lack  of  the  constraint-imposed 
shaping  effect,  the  control  histories  differ  significantly 
from  those  of  the  maneuver  enforcing  the  terminal 
orthogonality  constraint.  Mote  that  the  thrust  control 
angle  is  confined  to  thrusting  nearly  in  the  radial 
direction  over  the  entire  maneuver  for  the  two  hour  and  four 
hour  times  of  flight.  During  the  six  hour  time  of  flight, 
however,  the  thrust  angle  is  directed  into  the  second 
quadrant,  towards  a  terminal  position  again  on  the  +y-axis. 
Additionally,  the  two  and  four  hour  time  of  flight  maneuvers 
again  show  an  increased  orbit  period  at  their  termination, 
while  orbit  drift  is  eliminated  over  the  six  hour  time  of 
flight  maneuver.  The  effect  of  the  terminal  position  upon 
the  radial  increases  in  distance  obtained  is  similar  to  that 
of  the  maneuver  with  the  terminal  orthogonality  constraint, 
but  the  target  spacecraft  travels  a  much  greater  radial 
distance,  as  Table  VII  shows. 


Table  VII.  Radial  Increase  for  Maximum  Radius 
Evasive  Maneuver  with  no  Terminal  Constraint 


Time  of  Flight  (hrs)  Radial  Increase  (km) 

2  26.5151 

4  116.4321 


6 


259.0089 


TIME  (HOURS) 


Figure  9.  Thru9t  Control  Angle  V9  Time  for  2  Hour  Time  of 
Flight  to  Maximum  Radius  with  no  Terminal  Constraint 


TIME  (HOURS) 


Figure  10.  Thrust  Control  Angle  vs  Time  for  4  Hour  Time  of 
Flight  to  Maximum  Radius  with  no  Terminal  Constraint 


THRUST  CONTROL  ANGLE  (DEGREES) 


120.0 


Figure  1 1 .  Thrust  Control  Angle  vs  Time  for  6  Hour  Time  of 
Right  to  Maximum  Radius  with  no  Terminal  Constraint 


.ximum  radius  aantuvar .  Mot*  the  orientation  of  the 


apacacraft  velocity  vector  in  the  radially  outward 
direction. 


Table  VII I.  Terminal  Position  States  for  Maximum  Radius 
Evasive  Maneuver  with  no  Terminal  Constraint 


Time  of 

Flight 

(hrs) 

x  (km) 

y  (km) 

2 

36605. 1947 

21133.2181 

4 

21194.0656 

36673.8873 

6 

o 

• 

o 

42500. 1314 

Table  IX.  Terminal  Velocity  States  for 
Evasive  Maneuver  with  no  Terminal 

Maximum  Radius 
Constraint 

Time  of 

Flight 

(hrs) 

u  (km/s) 

v  (km/s) 

2 

-1.5292 

2.6640 

4 

-2.6488 

1.5513 

0 

-3.0710 

0.0323 

The  corresponding  control  histories  for  the  two,  four 
and  six  hour  time  of  flight  maneuvers  to  a  minimum  radius 
with  no  constraint  at  the  final  time  are  illustrated  in 
Figures  12,  13  and  14.  Again,  for  each  time  of  flight,  the 
thrust  angle  profile  is  nearly  identical  to  its  maximum 
radius  counterpart,  but  shifted  by  180  degrees.  The 
resulting  changes  in  radial  distance  are  then  only  in  a 
radially  opposite  direction,  as  shown  in  Table  Z. 


TIME  (HOURS) 


Figure  14.  Thrust  Control  Angle  vs  Time  for  6  Hour  Time  of 
Right  to  Minimum  Radius  with  no  Terminal  Constraint 


Table  X.  Radial  Decrease  for  Minima  Radius 
Bvasive  Maneuver  with  no  Terminal  Constraint 


Time  of  Flight  (hrs) 

Radial  Decrease  (km) 

2 

20.4681 

4 

116.4797 

6 

250.2153 

The  final  position  and  velocity  states  for  this  evasive 
maneuver  are  given  in  Tables  XI  and  XII. 


Table  XI.  Terminal  Position  States  for  Minimum  Radius 
Evasive  Maneuver  with  no  Terminal  Constraint 


Tim  of  Flight  (hrs) 

x  (km) 

y  (km) 

2 

30558.0342 

21107.8008 

4 

21047.1700 

30480.7502 

0 

O 

• 

o 

41081.0071 

Table  XII.  Terminal  Velocity  States  for  Minimum  Radius 
Evasive  Maneuver  with  no  Terminal  Constraint 


Time  of  Flight  (hrs)  u  (km/s)  v  (km /a) 


2  -1.5427  2.0560 

4  -2.0716  1.5205 

0  -3.0728  0.0030 

Rendezvous  Maneuver.  Thrust  vector  control  histories 
were  generated  for  minimum  time  rendezvous  maneuvers  with  an 
initial  position  at  the  terminal  point  of  the  maximum  radius 


evasive  maneuvers  with  the  orthogonality  constraint. 


Figures  15,  16  and  17  show  the  control  angle  profiles  for 
the  rendezvous  from  the  two  hour,  four  hour  and  six  hour 


termination  points,  respectively.  For  all  rendezvous 
maneuvers,  the  thrust  control  angle  is  nearly  constant  for 
long  periods  of  time,  with  the  control  switching  direction 
by  180  degrees  at  approximately  the  sixty  percent  time  point 
of  the  maneuver.  Although  the  control  angle  is  not  bounded 
in  the  problem  formulation,  the  numerical  solution  yields  a 
control  angle  region  bounded  approximately  by  -30  and  160 
degrees  for  the  rendezvous  from  the  two  hour  evasive 
■aneuver ,  0  and  180  degrees  for  the  return  from  the  four 
hour  eaneuver ,  and  -SO  and  00  degrees  for  the  rendezvous 
from  the  six  hour  maneuver.  The  penalty  for  achieving 
increased  radial  distance  at  the  expense  of  lagging  behind 
the  nominal  orbit  is  evidenced  in  the  time  required  for 
rendezvous  from  the  two  hour  and  four  hour  time  of  flight 
terminal  points,  as  Table  XIII  shows. 

Table  XIII.  Time  Required  to  Rendezvous  from  Evasive 
Maneuver  with  Terminal  Constraint 

Time  of  Flight  for  Rendezvous  Time 

Evasive  Maneuver  (hrs)  (hrs) 

2  3.6104 

4  8.0241 


5.0057 


TIME  (HOURS) 


Figure  15.  Thrust  Control  Angle  vs  Time  for  Rendezvous  from  2  Hr 
Time  of  Flight  to  Maximum  Radius  with  Orthogonality  Constraint 


To  verify  that  the  terminal  states  of  the  rendezvous 


maneuver  do  indeed  match  those  of  the  nominal  orbit,  the 
point  of  rendezvous  is  given  in  Table  XIV.  Table  XV  shows 
the  terminal  velocity  states  at  the  rendezvous  point. 


Table  XIV.  Terminal  Position  States  for  Bendezvous  from 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of 
Evasive 

Flight  for 
Maneuver  (hrs) 

x  (km) 

y  (km) 

2 

4301 . 1276 

42021.5792 

4 

-42238. 1810 

-261.9727 

6 

-42241 . 1149 

47. 1797 

Table  XV.  Terminal  Velocity  States  for  Rendezvous  from 
Evasive  Maneuver  with  Terminal  Constraint 


Time  of 
Evasive 

Flight  for 
Maneuver  (hrs) 

u  (k m/s) 

v  (km /s) 

2 

-3.0559 

0.3128 

4 

0.0191 

-3.0716 

6 

-0.0034 

-3.0719 

The  total  amount  of  time  required  for  a  complete 
evasion  and  rendezvous  maneuver  is  given  in  Table  XVI.  Here 
the  penalty  of  lagging  behind  the  nominal  orbit  in  order  to 
achieve  a  larger  radial  change  is  seen  dramatically.  Mote 
that  while  the  maneuvers  associated  with  the  4  hour  time  of 
flight  and  the  six  hour  time  of  flight  evasive  maneuvers 
require  approximately  the  same  total  flight  time,  a 


trade-off  exists  between  achieving  a  greater  relative  radius 


change  during  the  evasive  maneuver  and  returning  more  quickly  to 
the  nominal  orbit. 


Table  XVI.  Total  Time  Required  for  Evasive  and 
Rendezvous  Maneuver 


Time  of  Flight  for 

Evasive  Maneuver  (Hours) 

Maneuver  Time 
(Hours) 

2 

5.6104 

4 

12.0241 

6 

11.9057 

Rendezvous  following  any  minimum  radius  evasive 
maneuver  was  not  a  concern,  since  the  interceptor  was 
assumed  to  be  in  a  direct  coasting  ascent.  The  minimum 
radius  trajectory  could  possibly  pass  within  striking 
distance  of  the  ascending  threat  depending  upon  the  size  of 
its  lethal  radius,  which  was  not  specified.  However,  it  is 
anticipated  that  a  return  to  the  nominal  orbit  from  a 
minimum  orbit  radius  would  require  less  time  since  the 
trajectory  would  traverse  the  interior  of  the  nominal  orbit. 

On  the  other  hand,  rendezvous  from  a  maximum  radius 
maneuver  with  no  terminal  constraint  is  of  great  concern, 
since  this  would  make  apparent  the  penalty  paid  in 
rendezvous  time  in  returning  from  a  much  larger  radial 
distance  with  an  unconstrained  orientation  of  the  spacecraft 
velocity  vector.  Unfortunately,  convergence  could  not  be 


achieved  within  the  limits  of  the  8  megabytes  of  computer 
memory  available;  these  maneuvers  demanded  grid  sizes  in 
excess  of  2342,  the  upper  limit  available  at  this  time  on 
the  Vax- 11/785  system,  as  indicated  by  the  DVCPB  parameter 
MGBID. 

The  difficulty  posed  herein  lies  not  only  in  the  great 
return  distance  demanded,  but  also  in  an  initial  orientation 
of  the  spacecraft  velocity  vector  in  the  outward  radial 
direction.  In  fact,  convergence  was  achieved  for  rendezvous 
initiating  at  the  desired  x-y  position,  but  with  the 
velocity  vector  orthogonal  to  the  radius  vector.  This 
solution  was  intended  as  an  input  grid  for  the  return 
maneuver  from  the  unconstrained  evasive  maneuver  with 
non-orthogonal  velocity  and  radius  vectors.  In  order  to 
achieve  the  desired  solution,  the  initial  velocity  vector 
was  then  rotated  to  its  non-orthogonal  orientation  using  a 
continuation  scheme  with  small  step  size.  At  most  the 
vector  could  be  rotated  sixty-percent  towards  the  desired 
direction  before  computer  memory  limits  were  exceeded,  and 
even  then  it  was  apparent  that  the  rendezvous  time  required 
was  great  in  comparison  to  those  cited  in  Table  V. 

Minimum  Time  Intercept  Maneuver .  The  maximum  radius 
evasive  maneuvers  were  also  considered  in  the  context  of 
a  minimum  time  intercept  of  a  specified  radius,  with  and 
without  a  terminal  orthogonality  constraint.  The  resulting 


state  and  control  histories  duplicated  those  of  the  maximum 
radius  maneuvers,  thus  verifying  the  solutions. 

Costate  Trajectories 

As  discussed  in  Chapter  VI,  the  costate  trajectories 
vary  greatly  depending  upon  the  initial  grid  input  to  the 
problem  solver  DVCPB.  Figures  18  through  47  evidence  this 
behavior,  as  the  x-y  and  u-v  costate  trajectories  show  only 
limited  similarity  in  their  time  histories,  even  for  the 
same  maneuver  over  a  different  time  of  flight.  The 
exceptions  herein  are  the  costate  trajectories  associated 
with  mini mum  orbit  radius  maneuvers,  shown  in  Figures  24 
through  20  and  36  through  41  for  constrained  and 
unconstrained  maneuvers,  respectively.  These  trajectories 
are  merely  reversed  in  sign  from  their  max i mum  radius 
counterparts,  as  the  problem  development  suggests. 
Additionally,  the  u-v  costate  histories  for  all 
maneuvers  without  constraint  clearly  satisfy  the  boundary 
conditions  requiring  these  costates  to  pass  through  zero  at 
the  final  time.  For  all  maneuvers,  the  costate  trajectories 
are  smooth  and  continuous  over  time,  as  required  by 
Euler-Lagrange  theory. 

Although  the  costate  trajectories  show  only  limited 
trends  towards  similar  histories,  they  can  be  fit  quite 
accurately  with  appropriate  quadratics.  Tables  XVII  through 
XIX  give  the  quadratic  curvefits  for  x-y  and  u-v  costate 


trajectories  for  each  maximum  radius  and  rendezvous 


f 


maneuver,  along  with  the  corresponding  correlation 
coefficients.  Mote  that  curvefits  on  the  coefficients  for 
minimum  radius  maneuvers  are  merely  the  negative  of  the 
maximum  radius  curvefits.  The  purpose  in  listing  these 
equations  and  their  associated  plots  is  to  aid  the  future 
user  in  constructing  an  input  grid  which  will  readily  alio 
convergence.  In  each  equation,  the  time  considered  is  in 


hours . 


TIME  (HOURS) 


Figure  20.  X  and  Y  Costates  for  4  Hour  Time  of  Flight  to 
Maximum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  21.  U  and  V  Costates  for  4  Hour  Time  of  Flight  to 
Maximum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  23.  U  and  V  Costates  for  6  Hour  Time  of  Flight  to 
Maximum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  25.  U  and  V  Costates  for  2  Hour  Time  of  Right  to 
Mimimum  Radius  with  Orthogonality  Constraint  at  Rnal  Time 


Figure  26.  X  and  Y  Costates  for  4  Hour  Time  of  P:ght  to 
Minimum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  27.  U  and  V  Costates  for  4  Hour  Time  of  Right  to 
Minimum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  28.  X  and  Y  Costates  for  6  Hour  Time  of  Right  to 
Minimum  Radius  with  Orthogonality  Constraint  at  Final  Time 
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Figure  31.  U  and  V  Costates  for  2  Hour  Time  of  Flight  to 
Maximum  Radius  with  no  Terminal  Constraint 


Figure  32.  X  and  Y  Costates  for  4  Hour  Time  of  Right  to 
Maximum  Radius  with  no  Terminal  Constraint 


Figure  33.  U  and  V  Costates  for  4  Hour  Time  of  Right  to 
Maximum  Radius  with  no  Terminal  Constraint 


Figure  34.  X  and  Y  Costates  for  6  Hour  Time  of  Right  to 
Maximum  Radius  with  no  Terminal  Constraint 
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Figure  35.  U  and  V  Costates  for  6  Hour  Time  of  Right  to 
Maximum  Radius  with  no  Terminal  Constraint 
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Rgure  36.  X  and  Y  Costates  for  2  Hour  Time  of  Right  to 
Minimum  Radius  with  no  Terminal  Constraint 
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Figure  38.  X  and  Y  Costates  for  4  Hour  Time  of  Flight  to 
Minimum  Radius  with  no  Terminal  Constraint 
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Figure  40.  X  and  Y  Costates  for  6  Hour  Time  of  Right  to 
Minimum  Radius  with  no  Terminal  Constraint 
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Rgure  41.  U  and  V  Costates  for  6  Hour  Time  of  Right  to 
Minimum  Radius  with  no  Terminal  Constraint 
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Figure  43.  U  and  V  Costates  for  Rendezvous  from  2  Hour  Time 
of  Flight  to  Maximum  Radius  with  Orthogonolity  Constraint 


^  Y  Costate 


V  Costate 


Y  Coslate 


X  Costate 


TOF  (Hr a) 


Coatata 

Equation 

Correlation 

X 

X 

1 . 022- . 1590t+ . 0687 t* 

0.9996 

X 

y 

.3149- . 0305 t+ . 01 16t2 

0.8279 

5.210-4.318t+.0502tZ 

1.0000 

m 

.7880-1 .379t+.0235t2 

1.0000 

X 

X 

1 . 157- . 3487t+ . 0689t2 

0.9971 

X 

y 

. 7275- . 0858t+ . 0 143t* 

0 . 7674 

X 

u 

11. 11-4. 186t+. 1610t2 

0 . 9989 

X 

V 

5.635-3. 154t^.0897t2 

0.9999 

X 

X 

- . 0595- . 0540t+ . 0149t2 

0.9132 

X 

y 

. 4957+ . 035 1 t+ . 0007 t* 

0.7674 

X 

u 

- 1 . 221-  .  4993 1-*- .  0240 1* 

0.9904 

X 

V 

5 . 798-2 . 294t- . 0520t2 

0.9999 

Table  XVIII.  Coatate  Quadratic  Fits  for  Maximum  Radius 
Maneuver  with  no  Terminal  Constraint 


TOF  (Hr s) 


Costate 


Equation  Correlation 

1. 135-.2717t+.0685t2  1.0000 
. 5067- . 0290t+ . 0182t2  0.8308 
8 . 483-4 . 828 t+0 . 304 12  1.0000 
4 . 568-2 . 298 t+0 . 000 t2  1.0000 


1 . 518- . 5288t+ . 0685t2  0.0097 
.9171-. 1100t-.0346t2  0.8238 
12.83-3. 261t~. 1224t2  0.0715 
17 . 06-4 . 380t+0 . 000t*  0.0003 


.6474-.  lOBOt-.OHOt2  0.0884 

.4630+.4770t-.0748t2  0.0084 

-. 8232-3. 305t+.6S81t2  0.0078 

28. 19-3. 500t-. 2820t2  0.0086 


I  J.'tf  .  t<.  |»,  |*>  l  «.<  t.t  fc-»  >.<  ^  <»<.t  4-*  l.l‘l 


Tab la  ZIZ.  Costate  Quadratic  Fits  for  Mini bub  Til 

Rendezvous  Maneuver 


TOF  (Hrs)  Costate 


▼III.  Conclusions  and  Recommendations 


Orbital  avasiva  maneuvers  enabling  a  target  spacecraft 
in  geosynchronous  orbit  to  evade  an  interceptor  in  direct 
coasting  ascent  have  been  examined.  A  continuous  low  thrust 
propulsion  system  consisting  of  ten  30  centimeter  primary 
propulsion  thrusters  propelled  the  target  craft  over  each 
avoidance  maneuver.  The  thrust  vector  control  history  was 
determined  over  time  to  generate  the  optimal  spaceflight 
trajectory  achieving  a  maximum  change  in  orbit  radius.  In 
the  context  of  optimal  control  theory,  this  change  in  radius 
was  achieved  by  both  maximizing  and  minimizing  the  orbit 
radius  of  the  target  spacecraft. 

Two  strategies  were  considered  for  the  evasive 
maneuvers.  An  orthogonal! ty  constraint  was  first  applied  at 
the  final  time  of  each  maneuver  Specifying  the  target 
spacecraft’s  terminal  velocity  as  necessarily  orthogonal  to 
its  terminal  radius.  The  constraint  was  then  removed  such 
that  the  spacecraft  traveled  without  regard  for  the  final 
orientation  of  its  velocity  vector.  When  the  spacecraft  was 
allowed  to  travel  free  of  constraint,  the  terminal  thrust 
control  angle  became  indeterminate.  This  singularity  was 
resolved  by  applying  L’Hospital's  rule  to  the  bilinear 
tangent  control  law,  such  that  the  control  definition 
applied  to  the  state  equations  was  changed  at  the  final  time 
of  the  evasive  maneuver. 
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For  both  constraint  strategies,  optimal  saxlnua  radius 
and  optimal  mini mum  radius  trajectories  achieved  the  sane 
change  in  orbit  radius,  but  in  the  opposite  direction.  As  a 
result  of  comparing  the  two  classes  of  maneuvers,  it  was 
apparent  that  a  significant  penalty  was  paid  in  terms  of 
radial  distance  traveled  when  a  terminal  constraint  was 
applied.  A  second  comparison  was  available  in  that,  in  the 
case  of  either  strategy,  change  in  radius  was  sacrificed  for 
the  elimination  of  orbit  drift  to  retain  the  nominal  orbit 
period.  As  expected,  longer  time  of  flight  maneuvers 
yielded  greater  changes  in  orbit  radius.  It  must  be 
considered,  however,  that  the  time  of  flight  of  an 
implemented  maneuver  should  be  as  short  as  possible  to 
maintain  a  minimum  of  fuel  consumption. 

Minimum  time  intercepts  of  a  specified  radius  were 
considered  as  a  means  of  duplicating  the  maximum  radius 
evasive  maneuvers.  The  resulting  time  optimal  state  and 
control  trajectories  were  indeed  identical  to  results  found 
earlier,  thus  verifying  the  quality  of  the  maximum  radius 
solutions . 

Opon  completion  of  the  evasive  maneuver,  the  target 
spacecraft  was  required  to  rendezvous  with  its  nominal  orbit 
position.  Rendezvous  maneuvers  could  only  be  completed  for 
maximum  radius  maneuvers  with  the  terminal  orthogonality 
constraint  applied  due  to  computer  memory  limits.  In 
comparing  the  maneuvers  that  were  successfully  completed,  it 
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For  both  constraint  strategies,  opt i sal  saxisua  radius 


mas  evident  that  the  tine  required  to  rendezvous  was,  in 
general,  significantly  greater  than  the  tine  of  flight  of 
the  evasive  naneuver .  One  case  was  found,  however,  in  which 
a  six  hour  orbit  evasive  naneuver  enabled  an  approxinately 
equal  tine  of  flight  for  the  rendezvous  naneuver.  This 
anonaly  was  attributed  to  the  position  of  the  target 
spacecraft  at  the  end  of  the  six  hour  evasive  naneuver, 
which  was  on  the  sane  radius  vector  as  the  nominal  orbit  as 
a  result  of  the  elimination  of  orbit  drift.  Thus,  a 
trade-off  was  apparent  in  which  the  optimal  evasive  maneuver 
did  not  necessarily  transition  readily  to  the  optimal 
rendezvous  maneuver . 

Rendezvous  maneuvers  following  minimum  radius  evasive 
maneuvers  were  not  a  concern  since  the  evasive  trajectory 
may  tend  to  pass  through  a  lethal  radius  centered  on  an 
ascending  interceptor,  depending  upon  the  interceptor’s 
particular  characteristics .  It  was  anticipated,  however, 
that  these  maneuvers  would  consune  less  time  than  the  return 
maneuvers  from  maximum  radius  evasions  due  to  their  shorter 
path  length. 

The  optimal  rendezvous  trajectory  from  the  evasive 
maneuvers  with  no  terminal  constraint  could  not  be 
determined  by  numerical  methods  due  to  memory  limits  on  the 
computer  system  used.  It  was  anticipated,  however,  that  the 
return  from  a  position  with  an  unconstrained  velocity  vector 
would  require  signif icantly  more  tine  than  that  from  an 


orientation  determined  by  the  orthogonality  constraint. 

In  order  to  resolve  the  intercept-rendezvous  problem 
completely  in  terms  of  total  minimum  fuel  expended,  the 
numerical  solution  for  the  return  maneuver  from  an 
unconstrained  evasive  maneuver  must  be  determined  on  a 
computer  system  with  greater  memory  capacity.  Current  plans 
include  adapting  the  associated  two  point  boundary  value 
problem  solution  on  the  ASD  Cyber  system. 

The  development  and  solutions  herein  are  precursors  to 
modeling  a  target  spacecraft  in  three  dimensions.  Various 
dynamics  of  the  interceptor  may  also  be  considered,  such 
that  the  threat  is  itself  maneuvering.  The  formulation  of 
the  optimal  control  problem  would  then  maximize  the  distance 
between  the  threat  and  the  target  spacecraft  at  the  nominal 
intercept  time,  where  the  states  would  be  given  in  terms  of 
the  vector  from  interceptor  to  target  and  its  first  time 
derivative.  Such  a  description  of  the  two  point  boundary 
value  problem  would  determine  if  out  of  plane  maneuvers  are 
potentially  sore  optimal  than  those  in  plane.  These  studies 
could  yield  practical  results  pertinent  to  future  spacecraft 
mission  planning  and  the  increasingly  relevant  intercept 
avoidance  problem. 


SUBBOUT I HE  FCMI (N.T.X.XP) 

DOUBLE  PBECISIOM  T ,X(M) ,XP(H) .MU.B.B3 ,B5 ,B .D .TH.MO.MD , 
•  DSC 
IITEOEB  ■ 

COMMON  /A/  MU  /B/  TH.MO.MD 
B=DSQBT(X(1)**2+X(2)**2) 

B=TH/ (MO-MD*T) 

DSC=DSQBT (X(7) **2+X(8) **2) 

D*X(7) *X(  1) +  X(8) *X(2) 

B3>B»»3 
B5SB**5 
XP( 1) =X(3) 

XP(2) =X(4) 

XP  (3)  =-MU*X ( 1 )  /B3*B*X (7)  /DSC 
XP (4) --MU*X(2) /B3+B*X(8) /DSC 
XP (5) =MU» (X(7) /B3-3 . 0d+00*D«X ( 1 ) /B5) 

XP (6) -MU* (X (8) /B3-3 . 0d+00*D*X(2) /B5) 

XP (7) =-X(5) 

XP (8) =-X(6) 

BETUBB 

END 


SUB BOUT I HE  FCHJ(H.T.X.PD) 

DOUBLE  PBECISIOB  T ,X(M) .PD(M.H) .MU.B.B3 ,B5 ,B7 ,B ,D ,TH, 
*  MO.MD.E.F.Q.DSC 
INTEGEB  N 

COMMON  /A/  MU  /B/  TH.MO.MD 
B=DSQBT(X( 1) **2+X(2) **2) 

B=TH/ (MO-MD*T) 

B3=B**3 

B5=B**5 

B7=B**7 

DSC=DSQBT (X (7) »*2+X(8) *«2) 

D=X(7)*X(1)+X(8)*X(2) 

E=X(7) *X(2) +X(8) *X( 1) 

F*3 . 0d+00*X(7) *X( 1 ) +X(8) *X (2) 

G=3 . 0d+00*X(8) *X (2) +X(7) *X( 1 ) 


PD ( 1 . 1) sO . 0 
PD (1,2) -0.0 
PD(1 , 3) =1 . 00d+00 
PD (1,4) *0 . 0 
PD ( 1 . 5) =0 . 0 
PD ( 1 ,0) s0 . 0 
PD( 1 ,7) *0.0 
PD(1,8)«0.0 


¥ 


m 


PD (2, 1) *0.0 
PD (2. 2) *0.0 
PD (2 , 3) *0 . 0 
PD (2 ,4) *1 . 0d+00 
PD(2,5) *0.0 
PD(2,0) *0.0 
PD (2. 7) *0.0 
PD(2 ,8) *0.0 


PD (3 , 1 ) *MD« (- 1 . 0d+00/B3+3 . 0d+00«l ( 1 ) ••2/B5) 

PD (3 , 2) *3 . 00d+00*lfU*X  ( 1 ) *X (2) /BS 

PD (3 ,3) *0.0 

PD (3 , 4) *0 . 0 

PD (3, 3) *0.0 

PD (3 ,0) *0 . 0 


PD (7 ,7) *0 . 0 
PD(7.8)*0.0 


PD (8, 1) >0.0 
PD(8 ,2) *0.0 
PD (8 ,3) 30.0 
PD (8 , 4) =0 . 0 
PD (8, 5) *0.0 
PD (8 ,8) *- 1 . 0d*00 
PD (8 , 7) *0 . 0 
PD(8,6) *0.0 

RETUBM 

KMD 


SOB BO DTI ME  FCMB(M.XA.ZB.F) 

DOUBLE  PBECISIOH  ZA(H) ,XB(H) ,F(M) . KS.KM.RB 
IMTEOEB  > 

COMMON  /COMV/  KM.KS 
RB=DSQBT (XB ( 1 ) »»2+XB (2) «*2) 

F ( 1 ) *XA ( 1 ) -KM 
F (2) *XA(2) 

F (3) =XA (3) 

F(4)*XA(4)-KS 

F(5)  =XB(1)*XB(3)  ♦XB(2)*XB(4) 

F  (6)  =ZB  (5)  *XB  ( 1 )  -XB  ( 1)  ««2/BB-XB (7) #XB (3) 

F  (7)  -XB  (6)  *XB (2)  -XB  (2)  **2/BB-XB (8)  «XB  (4) 

F  (8)  -XB  (8)  *ZB  (1)  -ZB  (7)  «XB (2) 


Appendix  B:  Subroutine*  FCMI .  FCBJ  end  FCMB  for  lltxiw 
Bxdluj  Maneuver  with  no  Terwinal  Constraint 


SUBBOOTIME  FCMI (M.T.X.XP) 

DOUBLE  PBKCISIOM  T.Z(M) ,XP(M> .MU . H . B3 ,B5 . B . D ,TH . MO . MD . 
•  DSC 

INTEGER  M 

COMMDM  /A/  MU  /B/  TH.MO.MD 
B=TH/ (MD-MD*T) 

DSC-DSQRT(X(7) «»2*X(8) ••2) 

D*X(7)  *X(1)+X(8)*X(2) 

B3*B««3 
B5>Bh9 
XP  (1)  *X(3) 

XP (2) *X(4) 

IF  (DSC  .MB.  0 . 0d+00)  THEM 
XP (3) *-MU«X( 1) /R3*B»X(7> /DSC 
XP (4) *-MU«X(2) /B3*B*X(8) /DSC 
ELSE 

DSC=DSQBT (X(S) ««2*X(0) **2) 

XP (3) s~MU*X( 1) /R3*B«X(5) /DSC 
XP (4) * -MU»X ( 2 ) /B3+B*X(0) /DSC 
EMDIF 

XP(5)=MU*  (X(7)/B3-3.0d+00»D»X(l)/R5) 

XP(O) =  MU» (X (8) /R3-3 . 0d+00«D*X(2) /B5) 

XP (7) *-X(3) 


•  < 


XP (8) *-X(C) 

BETUBM 

EMD 


SUBBOUT I ME  FCMJ (M ,T . X.PD) 

DOUBLE  PBECISIOM  T.X(M) ,PD(M,M) ,MU ,R ,R3 ,R5 ,R7 ,B ,D ,TH , 
•  MO.MD.E.F.G.DSC 
INTEGER  M 

COMMON  /A/  MU  /B/  TH.MO.MD 
B*DSQRT(X(l)ft»2+X(2)*«2) 

B=TH/ (MO-MD»T) 

B3*R>«3 

B5=B««5 

B7=B*»7 

DSC=DSQBT (X (7) **2+X(B) ««2> 

D*X(7) »X( 1) +X(8) *X(2) 

E»X(7)«X(2)*X(8)«X(1) 

F*3 . 0d+00»X(7) *X(1)+X(8)*X(2) 

G-3.0d+00«X(8)«X(2) +X(7)«X(1) 

IF  (DSC  .ME.  0.0d+00)  THEM 


PD (1,1)' 
PD (1.2) 


>0.0 

>0.0 


m 


in 


PD ( 1 1 3) *  1 . 00d+00 
PD (1,4) *0 . 0 
PD(1 ,5) *0.0 
PD(1,8)*0.0 
PD(1,7)*0.0 
PD( 1 ,8) *0. 0 

PD(2,1)*0.0 
PD (2 ,2) *0 . 0 
PD (2 ,3) *0 . 0 
PD(2,4)*1 . 0d+00 
PD(2 ,5) *0.0 
PD (2 ,0) *0 .0 
PD(2.7)*0.0 
PD (2 ,8) *0 . 0 

PD (3 , 1 ) *M0* (- 1 . 0d*00/B3+3 . 0d+00*X ( 1 ) *»2/B9) 

PD (3 ,2) *3 . 00d+00*MU«X( 1)*X(2)/B5 

PD (3 ,3) *0.0 

PD(3,4)*0.0 

PD(3 ,5) *0.0 

PD(3,8) *0.0 

PD (3 .7) *B«X(6) ••2/DSC* *3 
PD (3 , B) *-B»X(7) «X(B) /DSC«»3 

PD (4 ,1) *PD (3 , 2) 

PD  (4 . 2)  -MD* (- 1 . 0d+00/B3+3 . 0d*00*X (2) »«2/B9) 

PD (4. 3) *0.0 
PD (4 ,4) *0.0 
PD(4 ,9) *0 . 0 
PD(4.0)*0.0 
PD (4 ,7) *PD(3 ,8) 

PD (4 ,8) *B«X(7) **2/DSC**3 

PD (5 , 1) *MD» (-3 . 0d*00»F/B9+ 19 . 0d+00»X ( 1 ) »«2»D/B7) 

PD (9 , 2) *MO» (-3 . 0d+00«E/B5+ 19 . 0d+00»X( 1 ) *X(2) «D/B7) 

PD (9 ,3) *0 . 0 

PD (9. 4) *0.0 

PD (9 , 9) *0 . 0 

PD (9 ,0) *0 . 0 

PD (9 , 7) *-PD (3,1) 

PD (9,8) *-PD(3 ,2) 

PD (0 , 1 ) *PD (5 , 2) 

PD (6 ,2) «MD* (-3 . 0d+00*Q/B9+ 19 . 0d*00*X (2) ••2*D/B7) 

PD(0 ,3) *0.0 

PD (8 ,4) *0.0 

PD (0,9) -0.0 

PD (0,8) *0.0 

PD (8 ,7)  *PD  (9 ,8) 

PD(0 ,8) *-PD(4 ,2) 


PD(7.1)*0.0 
PD(7 , 2) *0 . 0 
PD(7.3)*0.0 
PD (7 ,4) *0.0 
PDC7 , 9) *-l . 0d+00 
PD<7,8)-0.0 
PD(7 ,7) *0 . 0 
PD(7 , 8) *0 . 0 

PD (8, 1) *0.0 
PD (8 ,2) *0.0 
PD<8.3)«0.0 
PD (8,4) *0.0 
PD (8. 5) *0.0 
PD (8,8)*-l. Od+OO 
PD (8 ,7) *0 . 0 
PD (8 ,8) *0 . 0 

ELSE 

DSC*DSQBT (X (5) »«2+X(0) **2) 

PD( 1 , 1) *0.0 
PD (1,2) *0.0 
PD ( 1 , 3) *  1 . 00d*00 
PD( 1 ,4) *0.0 
PD (1.9) *0.0 
PD( 1 ,8) *0.0 
PD( 1,7) *0 . 0 
PD ( 1 ,8) *0 . 0 

PD (2 , 1) *0 . 0 
PD (2 ,2) *0 . 0 
PD(2,3) *0.0 
PD (2 ,4) *1 .Od+OO 
PD(2 ,9) *0.0 
PD (2 ,8) *0.0 
PD (2 ,7) *0 . 0 
PD(2 ,8) *0.0 

PD (3 , 1 ) -MU* (- 1 . 0d+00/B3+3 . 0d+00*X( 1 ) **2/B9) 

PD (3 . 2) *3 . 00d+00*MU»X ( 1 ) »X (2) /B5 

PD (3 , 3) *0 . 0 

PD (3 ,4) *0.0 

PD (3 ,7) *0 . 0 

PD(3 ,8) *0.0 

PD (3 ,9) *B«X(8) «»2/DSC»»3 
PD (3 ,0) *-B»X<9) «X(0) /DSC»»3 

PD (4 ,1) -PD (3 ,2) 

PD (4 ,2) *MU* (-1 . 0d+00/B3*3.0d*00ftX(2) «*2/B9) 
PD(4 ,3) *0 . 0 
PD (4 ,4) *0.0 


PD (4 ,7) »0.0 

PD (4 ,8) *0 .0 

PD  (4  ,S)  *  PD  (3,0 

PD (4 ,0) *B*X(9) ••2/DSC* *3 

PD  (5 , 1 )  =*M0*  (-3 . 0d+00»F/B5+ 15 . 0d+00*X  ( 1 )  ■*2*D/B7) 

PD (5 ,2) *M0* (-3 . 0d+00*K/B5+ 15 . 0d+00»I( 1 ) *X (2) *D/B7) 

PD (5, 3) *0.0 

PD (5 ,4) *0 . 0 

PD(9 ,5) *0.0 

PD(9.8)«0.0 

PD (9, 7) *-PD(3, 1) 

PD (9 ,8) *-PD(3,2) 

PD (6,1) *PD  (9 ,2) 

PD (8 ,2) *MJ» (-3. 0d+00«G/B9+19. 0d+00»X(2) »»2»D/B7) 

PD(0 ,3) *0.0 

PD (6 ,4) *0 . 0 

PD (6 a  9) *0 . 0 

PD (8 ,0) *0.0 

PD (8 ,7) *PD(9 ,8) 

PD (8 ,8)  *-PD  (4 , 2) 

PD(7.1)«0.0 
PD(7.2)*0.0 
PD (7,3) *0 . 0 
PD (7 ,4) *0 . 0 
PD (7, 9) *-i . 0d+00 
PD (7 ,8) *0 . 0 
PD(7 ,7) *0.0 
PD(7 ,8) *0 . 0 

PD (8, 1)*0.0 
PD(8 ,2) *0.0 
PD(8,3) *0.0 
PD (8 ,4) s0 . 0 
PD (8 , 9) *0 . 0 
PD (8, 8) *-l . Od+OO 
PD(8,7)«0.0 
PD (8 ,8) *0 . 0 

KWDIF 

bktubb 


8  UB  BOUT  in  FCMB (M.XA.XB.F) 

DOUBLK  PRECIS 101  XA(E) .ZB (■) ,F (■) , KS .KM.BB . AL 

IHTEQEB  ■ 

COIMOM  /COW/  KM.KS 
COMMOB  /COEST/  AL 


RB«DSQRT(ZB(1)*»2+ZB(2)»*2) 

F(1)-XA(1) -KM 
F (2) *ZA(2) 

F(3)*XA(3) 

F(4)*ZA(4) -KS 

F (9) *AL» (XB( 1) *ZB (3) +ZB (2) »ZB (4) )  ♦  ( 1 . Od+OO-AL) *ZB(7) 
F (0) =XB (5) »ZB( 1) -XB( 1) »«2/RB-AL» (ZB (7) »ZB(3> ) 

F (7) =ZB (0) »ZB (2) -ZB (2) ••2/BB-AL* (ZB (8) *ZB (4) ) 

F (8) =AL* (ZB (S) »ZB ( 1 ) -ZB (7) *ZB (2) ) ♦ ( 1 . Od+OO-AL) «ZB (8) 

RET 


I 


Appendix  C:  Subroutines  FCMI .  FCMJ  and  FCNB  for 
Rendezvous  Maneuver 


SUBROOT I ME  FCMI (M.T ,X.XP) 

DOUBLE  PRECISION  T.X(M) .XP(M) .MU. B.R3.R5.B.D. TH.MO.MD, 
•  DSC 
INTEGER  M 

COMMON  /ACOM/  MU  /BCOM/  TH.MO.MD 
R-DSQRT (X(l)«*2+X(2)**2) 

B=TH/ (M0-MD*T*X(9) ) 

DSC=DSQRT (X(7) *«2+X(8) »«2) 

D=X(7)«X(1)+X(8)*X(2) 

R3-R*»3 

R5*R»»5 

XP  (1)  SX(3)  *X(0) 

XP(2)*X(4)«X(9> 

XP(3)*(-MU«X(l)/R3-B»X(7)/DSC)eX(9) 

XP(4) = (~MU*X(2) /R3-B*X(8) /DSC) *X(9) 

XP (5) =MU» (X(7) /R3-3 . 0d+00»D»X( 1) /R5) »X(9> 

XP (6) =MU» (X (8) /R3-3 . 0d+00»D«X(2) /R5) «X<9) 

XP (7) *-X(5) *X(9) 

XP(8) *-X(6) *X(9) 

XP (9) *0 . 0d+00 

RETURN 

END 


SUBROUTINE  FCNJ (N,T,X,PD) 

DOUBLE  PRECISION  T,X(N) , PD(N,N) , MU , R , R3 , R5 , R7 , B , D , TH , 
i  MO.MD.E.F.G.DSC 
INTEGER  N 

COMMON  /ACOM/  MU  /BCOM/  TH.MO.MD 
R«DSQRT(X( 1) »»2+X(2) ««2) 

B=TH/ (M0-MD«T»X(9) ) 

DSC-DSQRT(X(7) **2+X(8) *«2) 

R3sR**3 

R5=R«*5 

R7*R««7 

D*X(7)  »X(  I)  +  X(8)  «X(2) 

E=X  (7)  *X(2)+X(8)»X(1) 

F=3 . 0d+00»X(7)  *X(  1)  +X(8)  *X(2) 

G=3 . 0d+00*X(8) *X(2)+X(7)*X(1) 


PD(l.l)-0.0 
PD( 1 .2) *0.0 
PD (1.3) *X (9) 
PD (1,4) m0 . 0 
PD ( 1 , S) *0 . 0 
PD( 1 ,8) >0.0 
PD(1 ,7) *0.0 


m 


m 


m 


w 


PD(1 ,8) *0.0 
PD(1,9)*X(3) 


PD (2, 1) =0.0 
PD(2.2)=0.0 
PD(2 ,3) *0.0 
PD(2 ,4) *X(9) 

PD(2 ,5) =0.0 
PD(2,0)*O.O 
PD (2 ,7) =0 . 0 
PD (2 , 8) *0 . 0 
PD(2 ,0) =1(4) 

PD (3,1) SMU* (-1 . 0d+00/R3+3 . 0d+00*X( 1 ) **2/B5) *X(9) 
PD(3.2)=3.00d+00*lfU«X(l)*X(2)/R5*X(9) 

PD (3 ,3) =0 . 0 
PD (3 ,4) *0.0 
PD (3 , 5) =0.0 
PD (3, 6) -0.0 

PD (3 ,7) =-B»X(8) **2/DSC**3«X(9) 

PD (3 ,8) =B»X(8) *X(7) /DSC«*3«X(9) 

PD (3,0) =  -lHJ*X(  1) /B3-B*X(7) /DSC 

PD (4 • 1) *PD(3 , 2) 

PD (4 . 2) =HD» (- 1 . 04+00/B3+3 . 0d+00«X (2) ««2/R5) »X(9> 

PD(4 ,3) =0.0 
PD (4 ,4) =0 . 0 
PD (4 ,5) *0 . 0 
PD (4 ,8) =0.0 
PD(4.7)*PD<3.8) 

PD (4 ,8) *-B*X(7) *• 2/DSC* *3«X(0) 

PD (4 ,0) *-M0»X(2) /B3-B*X(8) /DSC 

PD (3 , 1) =MD* (-3 . 0d+00*F/B5+ 15 . 0d+00«X ( 1 ) «*2»D/B7) *X  (9) 

PD (5 , 2) =MD* (-3 . 0d+00*E/R5+ 15 . 0d+00*X( 1) *X(2) *D/B7) *X(9) 

PD (5 ,3) *0 . 0 

PD(5 ,4) *0.0 

PD (5 .5) *0.0 

PD(3 ,0) >0.0 

PD (5,7) *-PD(3 ,1) 

PD (5 .8) *-PD(3 ,2) 

PD (5 ,0) *MD* (X(7) /B3-3 . 0d*00*X (1) *D/R5) 

PD (0,1) *PD (5,2) 

PD (8,2) *M0* (-3 . 0d+00*G/R5* 15 . 0d+00*X(2) ««2*D/R7) «X(9) 

PD (0 ,3) *0.0 

PD(0,4)*O.O 

PD (8, 5) *0.0 

PD(0,0)*0.0 

PD (0,7) *PD(5,8) 

PD (0,8) *-PD (4 ,2) 

PD(0 ,9)  ‘■MU*  (X(8)  /B3-3 . 0d+00*X(2)  *D/B5) 
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» 


PD (7 , 1) =0.0 
PD (7, 2) =0-0 
PD(7,3)*0.0 
PD (7 ,4) =0 . 0 
PD(7,5)=-X(9) 
PD (7 ,6) =0  -  0 
PDC7 ,7) =0 . 0 
PD(7,8) =0.0 
PD(7,9)=-X<5) 

PD (8 , 1 ) =0 . 0 
PD (8 , 2) =0 . 0 
PD (8 ,3) =0 . 0 
PD (8,4) =0.0 
PD (8 ,5) =0 . 0 
PD (8 ,6) =-X (9) 
PD(8,7) =0.0 
PD (8 ,8) =0 . 0 
PD (8 ,9) =-X (6) 

PD (9,1) =0 . 0 
PD (9 , 2) =0 . 0 
PD (9 ,3) =0 . 0 
PD(9,4) =0.0 
PD (9 ,5) =0 . 0 
PD (9 ,8) =0.0 
PD (9 ,7) =0 . 0 
PD(9 ,8) =0.0 
PD (9 ,9) =0 . 0 

BETUBH 

BHD 


SUBBODTIME  FCMB (H.XA.XB ,F) 

DOUBLE  PRECISIOH  XA(H) ,XB (H) ,F (H) , KM.KS ,BB3 ,XBP3 ,XBP4 , 

•  TH,M0,MD,MU,BF,XBM0M(4)  .DCB.W.TP.AL 
IHTEGEB  H 

COMMOM  /ACOM/  MU  /BOOM/  TH.MO.MD  /COHV/  KM.KS  /DCOM/  TP 
COMMOH  /COHST/  AL 
W=KS/KM 
BBS = KM* *3 

BF-TH/ (M0-MD«XB (9) ) 

DCB=DSQBT (XB (7) *«2*XB (8) **2) 

XBP3=XB (7) • (~MU«XB ( 1 ) /BB3-BF*XB (7) /DCB) 

XBP4=XB (8) • (~MU*XB (2) /BB3-BF«XB (8) /DCB) 

XBHOM( 1 ) =KM«DCOS (M«XB(9) +TP) 

XBH0M(2)  =KM»DSIM(W«XB(9)  +TP) 
XBM0M(3)=-KS«DSIM(W«XB(9)+TP) 

XBH0M(4) =KS»DC0S (W*XB(9) +TP) 
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F ( 1 ) =XA( 1 ) - . 332645d+0 1 
F(2) =XA(2) - . 573847d+01 
F (3) =XA (3) ♦ . 335233d+00 
F (4) =XA (4) - . 194326d+00 

F (5) = ( 1+XBC5) *XB (3) +XB(6) *ZB(4) +XBP3+XBP4) «XB(9) 

»  ♦ (XB(5) *XBH0M(2) -XB(6) *KBM0M(1) +XB(7) *ZBH0M(4) 

i  -KB (8) *XBV0M(3)  )  *W 

F (6) =XB ( 1 ) +4 . 54425d+00«AL+ ( AL- 1 . Od+OO) «XBMOM( 1 ) 

F (7) =XB (2) -4 . 81623d+00»AL+ (AL- 1 . Od+OO) •XBM0M(2) 

F (8) =XB (3) + .282583d+00»AL+ (AL-1 .Od+OO) *XBHOM(3) 

F (9) =XB (4) ♦ . 266625d+00*AL+ ( AL- 1 . Od+OO) «XBHOM(4) 

RETURN 

EHD 
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Optimal  in-plane  orbital  maneuvers  enabling  a  target  spacecraft  in 
geosynchronous  orbit  to  avoid  interception  by  a  threat  in  direct  coasting 
ascent  axe  examined.  Given  constant  magnitude,  low  thrust  electric 
propulsion,  the  thrust  vector  control  history  is  determined  to  propel 
the  spacecraft  to  either  a  maximum  or  minimum  orbit  radius  in  a  specified 
time  of  flight.  The  evasive  maneuver  is  studied  from  two  aspects. 

First,  an  orthogonality  constraint  is  applied  at  the  final  time  of  the 
maneuver,  such  that  the  target  spacecraft's  position  vector  is  orthogonal 
to  its  velocity  vector.  This  requirement  is  then  removed  to  allow 
the  spacecraft  to  travel  without  influence  of  a  terminal  constraint. 
Results  demonstrate  the  adverse  shaping  effect  the  orthogonality 
constraint  applies  over  the  entire  evasive  maneuver. 

Upon  completion  of  the  evasive  maneuver,  the  spacecraft  returns 
in  minimum  time  to  its  nominal  orbit,  in  a  position  as  if  no  evasive 
thrusting  had  occurred.  Results  show  the  time  required  to  rendezvous 
with  the  nominal  orbit  may  be  significantly  greater  than  the  time  of 
flight  of  the  evasive  maneuver. 

Optimal  control  is  applied  to  formulate  the  evasive  and  rendezvous 
maneuvers  using  Euler-Lagrange  theory  and  the  calculus  of  variations. 

The  behavior  of  the  Lagrange  costates  is  examined  in  terms  of  their 
effect  upon  solution  convergence.  The  resulting  two  point  boundary 
value  problem  is  solved  by  numerical  methods  to  yield  the  optimal 
spaceflight  trajectory.  ( ~rl  <  sc  i  )  ^ 
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